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^ ■ ABSTRACT 

o : 

■ Aims. The purpose of the present investigation is to probe the dust properties inside a molecular cloud, in particular how particle 

^ ■ growth and the presence of ice coatings may change the overall shape of the extinction curve. 

^ ' Methods. Field stars behind a molecular cloud can be used to probe the cloud extinction for both the reddening and the absorption 

Ql features. By combining multi-colour photometry and IR spectroscopy the spectral class of the star can be determined as can the 

extinction curve, including the vibrational bands of ices and silicates. 
y—i • Results. Based on observations of field stars behind the dark globule B 335, we determine the reddening curve from 0.35 to 24 fjm. 

The water ice band at 3.1 yum is weaker (r(3.1) = 0.4) than expected from the cloud extinction (Av » 10 for the sightline to the most 
I I ' obscured star). On the other hand, the CO ice band at 4.7 yum is strong (t(4.67) = 0.7) and indicates that the mass column density of 

, frozen CO is about the same as that of water ice. We fit the observations to model calculations and find that the thin ice coatings on 

Othe silicate and carbon grains (assumed to be spherical) lower the optical extinction by a few percent. We show that the reddening 
' curves for the two background stars, for which the silicate band has been measured, can be accurately modelled from the UV to 

(— I I 24 fim. These models only include graphite and silicate grains (plus thin ice mantles for the most obscured star), so there is no need 

for any additional major grain component to explain the slow decline of the reddening curve beyond the K band. As expected, the dust 
I . model for the dense part of the cloud has more large grains than for the outer regions. We propose that the well established shallow 

O ' reddening curve beyond the K band has two different explanations: larger graphite grains in dense regions and relatively small grains 

^ ' in the diffuse ISM, giving rise to substantially less extinction beyond the K band than previously thought. 

C/2 , Conclusions. For the sightline towards the most obscured star, we derive the relation A^^j = 0.97-E(J-Ka:j), and assuming that all 

. silicon is bound in silicates, N(2H2+H) a 1.5-10^'-Av a 9- 10"' -A^rj. For the rim of the cloud we get Aj:, = 0.51 ■E(J-Kj), which is 

close to recent determinations for the diffuse ISM. The corresponding gas column density is N(2H2-l-H) a 2.3-l(f -Ay ~ 3-10^^-As:j. 

^-H ' Key words. B335 - Bok globule - interstellar extinction 

> ■ 

' 1. Introduction structed a dust model that describes the observed properties 

„ , , . , , , . , , for the diffuse ISM. It basically consists of silicate and carbon 

^ Interstellar reddening is caused by small dust particles, but .^^^ molecular clouds, the lower temperature, the higher 
the detailed optical and chemical properties of these particles ^ ^^^^^^ ^.^^i^^i^j^ ^^j^ ^jj^^ -^^^ ^^^^ 

^ remain largely unknown in particular for molecular clouds. ^^^^^^ condense on the grains. In addition, the grains coagulate 
^ With the exception of the 220 nm extinction peak (probably ^.^^ distribution is changed towards larger grains. These 

^ caused by small graphite grains), the UV, optical and near-IR ^Her the optical properties of the grains. The most 

> . spectral regions offer little information about the dust particles ^^^.^^^ difference compared to the diffuse ISM is the absorption 
except that they must be small compared to the wavelength ^^^^^^^^ ^^^^ ^^^^ ^j^^ instruments provided extensive 
(even in the UV). In the mid-IR, however, strong absorption ^roscopic evidence of various ice absorptio n features, 

^ . bands show up, the most universal of these being the broad band ^^^^^^^ h,0, CO, CO2, CH4, CH3OH and NH3 (lGibbetaI.1 

■ - . centred at g jjum due to silicates. Another widespread compo- 2004). In most cases, these absorption features have been ob- 

nent in the diffuse ISM (interstellar medium) PAHs (polycychc ^^^^^^ ^^^^^^^ luminous, deeply embedded protostars, and it is 

aromatic hydrocarbons) are seen m the mid-IR as emission (due ^^^^ circumstellar matter (influenced by dynami- 

to luminescence). In addition, the far-IR thermal emission from ^^^.^^.^^ processes involved in the formation of the star) 

the interstellar dust grains gives further constraints on their size ^^.^^ unprocessed cloud component. In addition, the intrinsic 

distribution, and to some extent, their composition. ^^j^ ^^p^^^^^j ^^^^.^^ distribution) of an embedded young stai- 

, . . , . . r ■ IT^ — : — I L , is not in general well known, which means that it is difficult to 

Combining all this info rmation, ^ramd (|2Q03b|) has con- accurately separate interstellar reddening from intrinsic IR ex- 

* Based on observations collected at the European Southern cess. For these reasons, field stars behind a molecular cloud are 

Observatory, Chile (ESO programs 073.C-0436(A), 077.C-0524(A)) better probes in investigations of the (more or less) undisturbed 

and on observations made with the Nordic Optical Telescope (NOT), cloud. The problem is to find background stars that are heavily 

operated on the island of La Palma jointly by Denmark, Finland, obscured to probe the dense regions and still bright enough for 

Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque IR spectroscopy. In spite of these difficulties, there are some 
de los Muchachos of the Instituto de Astrofisica de Canarias. 
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inv estigations of ice features using background stars: Elias 16 
by iGibb et all (l2004l) and Spitzer ice-data from three sources 
behind T aurus and CK2 i n Serp ens by Knez et al. (2005). More 
recently, iBoogert et al.l (1201 ll) has combined ground-based 
and Spitzer spectroscopy data to analyse the ice and silicate 
features in quiescent clouds by observing background field stars. 




Apart from the absorption bands, the continuous extinc- 
tion in the mid-IR for molecular clouds has been the s ubject 
of several recent investigations (e.g. 'Flaherty et al.' '20()7t 
Roman-Zunigaet al. 2007; Chapman et al. 2009; Boogertetal] 
201 II) . However, it is important to note that what is actually 
determined is the reddening curve, and to establish the extinc- 
tion curve, an assumption is made on the extinction ratio for 
two wavelengths, e.g. A h/Ak = 1.55 , which goes back to the 



frequently cited paper bv lRieke & Le bofskv ( 1985). Their paper 
combined earlier investigations with mid-IR observations of 
distant field stars of unknown intrinsic colours. Obviously the 
data available for their analysis was very limited compared to 
those in recent investigations, based on large surveys both from 
ground and space. The reddening curve is relatively straight 
forward to determine, either by statistical methods (slopes in 
colour/colour diagrams) or by classifying selected background 
stars and compare the observed colour to the intrinsic. By 
observing stars all the way to Spitzer/MIPS 24 fim one would 
expect that the extinction is so close to zero at this wavelength, 
that it for all practical purposes does not matter if it is assumed 
to be exactly zero. However, the observations show that the 
reddening curve levels off in the mid-IR, which is not expected 
from current dust models, and for this reason the extrapolation 

to zero wave numbe r remains uncertain. 

In our recent paper dOlofsson & Olofssonll20ld hereafter called 
Paper II) we confirmed that the reddening curve for molecular 
clouds in the wavelength range 0.35-2. 2 yum can be c harac- 
terized in a functional form described by ICardelli et al.l (Il989l) 
(henceforth CCM) with only one parameter, Ry = Av/Eb-v- It 
was, however, clear that the CCM function did not adequately 
describe the reddening curve beyond the K band, and in the 
present paper we extend the investigation into the mid-IR spec- 
tral region. For our case study we have compiled our own and 
archive data for a few stars behind the B 335 globule including 
photometry and spectroscopy from the optical region to 24 fj.m. 
Our main purpose is to characterize the mid-IR extinction. In 
addition we explore whether both the reddening curve, from 
UV to 24jL/m, and the absorption features can be explained in a 
single model. 




Fig. 1. A U band image of the B335 globule. The centre YSO and 
its outflow region as well as the background stars are marked. 
The circle represents the position of the '^CO an d the C^^O mea- 
surem ent closest to the background star #947 in lHariunpaa et al.l 
(I2OOI and the size of the circle represents the 45" beam width. 



Spectroscopy of the background star (# 947) in the water ice 
band at 2.7 < A < 3.8 //m and the CO-ice band at 4.0 < 
A < 4.5 yUOT using VLT-ISAAC at ESO, Paranal during two 
nights, 2004-04-12-13. A standard star HR7519 is observed 
before and after each observation of the program star. 
The 2D-spectra are first corrected for bias and flat-field 
and the ID-spectra are extracted with the IRAF program 
pipeline 'APALL'. The program and standard spectra are 
then divided by respective stellar atmospheric models. 
These are then divided to correct for telluric atmospheric 
extinction. 

Imaging at 3.79 fim using SIRCA (Stockholm Observatory 
IR CAmera) mounted on the 2.5 m Nordic telescope on La 
Palma. 



2. Observations and reductions 

The present investigation is based on observations of a few back- 
ground stars behind the B 335 dark globule (RA(2000) = [19 
37 00], DEC(2000) = [7 34 10]) in both the optical and the 
near-IR. Their positions are marked in Fig[T] Table |2] shows the 
observations available in our investigation. 



- Imaging and spectroscopy using the NTT/EMMI at ESO, 
La Sifla (077.C-0524(A) during four nights 2006-06-27-29. 
These obse rvations and reductions are d escribed in a previ- 
ous article, lOlofsson & OlofssonI ( |2009|) (Paper I). 



These data have been combined with data from the ISO and 
the Spitzer archives as summarized in Tables[T]and|2] 

3. Results 

3.1. The intrinsic colours of tlie baci^ground stars 

In Paper II we confirmed that the molecular cloud reddening in 
the optical and the near-IR spectral region is well characterized 
by a functional form defined by CCM with only one parameter, 
Rv = Av/Eb-v- This also allows us to use the multi-band pho- 
tometry to both determine the Ry value and the spectral class for 
an obscured star as well as the reddening towards it. For the two 
most obscured stars in the present study, #947 and #90, we get 
the same Ry value, 7?!/ - 4.9 (Fig. |2]i, and their spectral classes 
are M6III and AOV respectively. The spectral class of the most 
obscured star in our sample, #947, is of particular importance 
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Table 1. Observations and data. 



observation 


band 


observatory 


wavelengths jim 


observation date 


reference 


imaging 


optical 


NTT-EMMI 


[0.354, 0.413, 0.509, 0.673, 0.798] 


2006-06-27-29 




spectroscopy 


optical 


JN i 1-nJVlMl 


U.J < A < U.O 


onriA r\/^ t7 on 
zUUo-Ud-z /—Zy 




imaging 


NIR 


2MASS 


[1.258, 1.649,2.157] 




SknitskieetaL (2006) 


imaging 


L' 


NOT-SIRCA 


3.79 


2003-05-17 


unpublished 


spectrometry 


HjOice 


VLT-ISAAC 


2.55 <A< 4.2 


2004-04-12-13 




spectrometry 


CO ice 


VLT-ISAAC 


4.45 < ^ < 5.1 


2004-04-12-13 




spectrometry 


silicate 


ISOCAM-CVF 


5 <A< 16.5 


1997-04-25 


Rev-TDT-OSN: 526 003 74 


spectrometry 


silicate 


Spitzer-IRS 


5 < ^ < 21 


2005-10-11 


AOR: 10705664 VeluCores5 


CO isotopologues 


mm 


ESO-SEST 


"CO and C"*0 1-0 


1993 - 1998 


Harjunpaii et al (2004) 



Table 2. Summary of the background stars and available mea- 
surements. 



central star background stars_| 



star# 


YSO 


2 


10 


20 


84 


90 


112 


947 


ra 


19 37 00 


19 36 57 


19 36 54 


19 37 01 


19 36 56 


19 36 57 


19 36 57 


19 36 58 


dec 


07 34 10 


07 35 23 


07 33 49 


07 32 43 


07 33 44 


07 35 16 


07 33 30 


07 33 59 


Tgff 




3800 


4000 


9200 


6000 


9800 


6400 


3050 


nh nf nm ptr'v 


PPTltT*fll 
veil LI <tl 






flux [m 


Jy] 










vT ll T ^1^11.1^ Vim. ILlI 1 1 \ 
















NTT-EMMI 


















U602 


0.354 


0.004 


0.010 


0.074 


0.003 


0.005 


0.003 




Bb605 


0.413 


0.12 


0.22 


0.35 


0.021 


0.035 


0.019 




g772 


0.509 


0.71 


0.92 


0.85 


0.08 


0.087 


0.06 




r773 


0.673 


11.5 


7.7 


4.8 


0.69 


0.54 


0.49 


0.026 


1610 


0.798 


23 


11 


6.1 


1.0 


0.75 


0.73 


0.34 


2MASS 


















J 


1.258 


257 


55 


29 


5.9 


4.8 


4.6 


73 


H 


1.650 


612 


85 


39 


10 


6.8 


7.2 


357 


K, 


2.157 


661 


78 


35 


9.3 


7.1 


6.5 


627 


SIRCA 


3.79 














453 


Spitzer 


















IRACl 


3.550 




38 


15 


4.3 


3.6 


3.0 




IRAC2 


4.493 


220 


22 


9.8 


2.8 


2.5 


2.0 




IRAC3 


5.731 


169 


16 


6.7 


1.5 


1.9 


1.3 


243 


IRAC4 


7.872 


93 


8.9 


3.5 


0.8 


1.2 


0.7 


149 


MIPS 


23.7 


9.9 


0.85 


0.2: 




0.1: 


0.1: 


35 


spectrometry 


wavelength 






spectral resolution R 










range \pm] 
















NTT-EMMI 


0.4 < /I < 0.8 


+ 


+ 


-1- 


R^400 


-1- 


-1- 


+ 


VLT-ISAAC 


2.55 <A<4.2 














R ^600 


VLT-ISAAC 


4.45 < ^ < 5.1 














R ^3300 


ISOCAM-CVF 


5 <A< 16.5 














100 < < 240 


Spitzer-IRS 


7 < .!< 14 


60<R< 120 














Spitzer-IRS 


14 < ^ < 21 


60<R< 120 












60 <R< 120 



" numbers refer to star positions in Fig[T] 



as it is the only star for which we have water and CO ice data 
and the temperature determined by the multiband photometry is 
confirmed both by the optical spectrum (Fig. |3]l and the L win- 
dow spectrum (Fig. |4]i. The spectral classes of the other stars 
are also confirmed by optical spectra. This means that we can 
characterize the intrinsic SEDs (spectral energy distributions) of 
the stars, and for that purpose we use model atmospheres from 
iHauschildt et"an ([l999). 

3.2. The H2O and CO ice bands 

By dividing the ISAAC spectra of star #947 (corrected for 
telluric absorption and relative flux calibrated) with a model 



atmosphere corresponding to the spectral class of the star (Teff 
- 3050 K and log g = 0) we determine the absorption of the 
H2O and CO ice bands (see Figs. |4] and |5] respectively). The 
peak position of the water ice band is at 3280 cm"' (= 3.05 fim) 
with an FWHM of 350 cm We note that there is no clear 
extension shortward of 3000cm"' (longward of 3.33 fim). Such 
an extension is often observed and explained as due to methanol 
ice (iGibb et al.l2004 . The CO absorption peaks at 2140cm"' (= 
4.67 A fm), which differs from that of pure CO ice (~ 2142 cm"', 
lElsila et al. 1997), but agrees with observations of other regions 
(iPontoppidan et al.l l2003h . There is, however, no indication of 
wing components, although our S/N ratio does not allow any 
firm statement on this point. 
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n n l I I I I I I I I I I I I I I I I I I I I I I I I 

0.5 1.0 1.5 2.0 2.5 

wavenumber[1 //im] 

Fig. 2. The optical and near-IR extinctions towards the back- 
ground stars #947 and #90. The extinction for #90 and #947 
are well represented by the fitted CCM-extinction with Ry - 4.9 
(dashed). 




6000. 7000. 8000. 9000. 

wovelength [A ] 



Fig. 3. The optical spectrum of the most obscured star (#947) is 
well represented by a stellar model (from lHauschildt et al.ll999h 
with Teff 3050 K. The stellar model (red curve) has been modi- 
fied by a weighting function oc /("'to account for the interstellar 
reddening 




0.28 0.30 0.32 0.34 

wavenurnber [1/;irTi] 



Fig. 4. The ISAAC spectrum in the L band of star #947 (blue 
curve) and the stellar atmospheric model (red curve) are divided 
to get the water ice band (black curve). 



0.6 - 




J I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.2130 0.2135 0.2140 0.214b 0.2150 0.2155 
waveNumber[1 / fim] 



Fig. 5. The ISAAC spectrum in the M band of star #947 (blue 
curve) and the corresponding stellar atmospheric model (red 
curve) are divided to get the CO ice band (black curve). 



3.3. The silicate absorption band 

The B 335 globule was observed with the ISOCAM/CVF and 
the background M giant (#947) is the only star visible in the 
field. The LW part of the scan consists of two spectral scans 
(back and forth) and we have adopted the average of the two. 
The wavelength coverage is 5 < /I < 16 /urn, and except for the 
longest wavelengths, the S/N is good, see Fig|6] 

Two of the background stars (#2 and #947) have also been 
observed by Spitzer/IRS in one partial band each. Fig |6] The 
Spitzer/IRS and the ISOCAM/CVF spectra of star #947 have 
been combined to one spectrum. We note that the band peaks 
at 0.109 yum-i (= 9.2 fiin) with an FWHM of ~ 380 /Jin. The 
silicate feature peak is thus found at a sho rter wavelengths 
than e . g. described by iDraine & Led d 19841) and iGibb et all 
(l2004l) . while the FWHM-widths agree with those found for star 
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0.25 - 



0J5 - 



0.05 - 




wavenumber [1/^m] 

Fig. 6. Spectra for sight lines towards stars #2 (lower) and #947 
(upper curve) as observed by Spitzer-IRS (#2: 7 < A < 14 jim) 
and ISOCAM-CVF (#947: 5 < A < 16 ^im) and photometry 
from Spitzer-IRAC and Spitzer-MIPS (+ signs). 



forming regions bv lGibbet al1(l2004l) . In a recent investigation, 
IChiar et al. (2011), who investigate absorption bands in the qui- 
escent cloud ICS 146, also show that the silicate band is slightly 
shifted towards a shorter wavelength. It is also interesting to 
note that even though #947 is much more obscured than #2, 
the silicate bands have approximate l y the same strength. This 
agrees with the results o f IChiaret all (1201 Ih . 



0.4 - 



■'i 0.2 - 



0.0 - 




wavenumber [ 1 /^m 



Fig. 7. The optical depth of the silicate band TsUicate towards star 
#2 (dashed line) and #947 (full line) as determined from Spitzer- 
IRS and ISOCAM-CVF observations, respectively. Although the 
S/N of the observations not allows for an interpretation of the 
details, we note that for star #947 the band peaks at w 0.109 
/im"' corresponding to 9.2 yum. 

that the scatter in the diagram reflects real qualitative variations 
of the extinction in the intervening cloud. Even though we can- 
not either exclude such variations, it is in our view meaningful to 
calculate an error- weighted average for the "rim" stars compared 
to that of the most obscured stars. The result is shown in Fig. |9] 
and we find that the extinction is only going down marginally 
between the Kj band and the silicate band. There is an apparent 
difference between the extinction in the "rim" region compared 
to the deeper part of the globule (represented by the direction 
towards star #947). In view of the uncertainties involved, it is 
not quite clear that this difference is real, but in the following we 
assume that it in fact is real and model the two cases separately. 



3.4. The CO2 band at 15.2 /im 

There is an indication in the CVF spectrum of the CO2 band 
at 15.2 fj.m but it is not confirmed in the IRS spectrum, which 
has a higher S/N ratio. We can therefore only give an upper 
limit Tcoi < 0.15, corresponding to a C02-column density 
Nco2 < 3 ■ 10'^ [cm-2]. 



3.5. The reddening curve in the mid infrared 

For each star in our sample we compare the observed colour 
indices to that calculated without interstellar reddening. We refer 
the colour excess to the K j band 

E(A - Ks) = (rriA - Ks)ohserved - (niA - Ks)inirin.sic 

and normalize by dividing with E(I - Kj). The result is shown 
in Fig. [8] where we also include the Spitzer spectra. The gen- 
eral impression is that there is not much reddening beyond the 
band, but in view of the uncertainties involved, both related 
to the observations and the model atmospheres, we cannot claim 



4. Interpretation 

4.1. Extinction and grain properties 
4.1 .1 . Grain size distribution 

Ma this et al] d 19771) (henceforth MRN) constructed grain 
models for the ISM extinction with two simple grain size 
distributions dN(a)/da for silicates and graphite respectively, 
both with dN(a)/da oc a"^-^ (with a as the grain radius) and 
0.001 < a < 0.25 fim. Several attempts to find the grain 
size distribution of the ISM from t he measured extinction have 
been done. Weingartner & Draind (12001.) have elaborated on 
this and suggested size distributions for graphite and silicate 
grains based on more than ten parameters and explaining a 
large wavelength range of the ISM-extinction. We followed 
their approach to analyze the B33 5 extinction curve in the op - 
tical and near-IR wavelength range dOlofsson & Olofssonl2010l) . 

The same approach in the wavelengths beyond 2 fim for 
a molecular cloud does not seem to be that successful. 
We have therefore taken a somewhat simpler approach as- 
suming a modified MRN grain size distribution still with 
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0.1 



0.0 



-0.1 



-TJ -0.2 



-0.3 



-0.4 




9471 



-ui 

CL 



o o o o 
< < < < 
K a: [T a: 



0.1 

wavenumber[1 //^n] 

Fig. 8. The "reddening" curve in the mid-IR, "normalized ex- 
cess" = E(A - Ks)/E{l - Ks) vs wave number, for six stars be- 
hind the B 335 globule. The normalized excess derived from the 
Spitzer/ISOCAM spectra of two of the stars, #2 and # 947, ai-e 
included with the same normalization. 



spherical graphite 
optica l refractive 



and silicate grains and st ill using the 
indices as proposed by iDraine & Led 



( 11984 and Draini ( l2003d) (in tabulated f ormat foun d at 
http://www.astro.princeton.edu/~draine/dust/du st.diel.html) l. 
We also increased the grain size range to 0.001 < a < 5 fim. 
The distributions are assumed to follow a power law grain size 
distribution dN/da (with the spherical grain radius a) for each 
substance (graphites and silicates): 

^^1 -r -P 
—;—\MRN - ^MRN ■ a 
da 

without any restrictions to the parameters. We have furthermore 
introduced an accumulation of grains with sizes around one 
grain-size by adding a gaussian grain distribution to each of the 
graphite and the silicate MRN-distributions. 



dN_ 

dN_ 
da 



\sauss ~ 1 . + Cp. 



dN. 



exp{ 



—(a — aY 



) 



width 



dN. 



, \MRN ■ , \eauss 

da da 



We use a x optimisation scheme to calculate the parameter 
sets Cmrn, Cgauss, d, aH-idih (for graphites and silicates 
respectively). We find that the two different extinction curves 
determined above actually can be accurately modelled by this 
kind of grain size distribution (Figs. [TOl and fTTI). Thus, there is 
no need to include any other major dust component than silicate 
and graphite grains. The corresponding size distribution of the 
grains are shown in Fig. [12] Even though we cannot draw too 
far-reaching conclusion about the details of the size distribu- 
tions, we note that both silicate and graphite grains tend to pile 
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Fig. 9. The normalized excess, (E(/l - Ks)fE,{\ - Kj), in the mid- 
IR for the weighted average of the five stars at the rim of the 
B 335 globule (blue ) and the most o bscured star (red). The green 
curve represents the lDrain3 (l2003al) Ry - 5.5 curve. 



up at a single grain size. However, the peak for the graphite 
grains is located around 2 fim for the more obscured star (#947) 
and at 0.5 //m for the star closer to the rim of the globule (#2). 
There is also a striking difference in the graphite/silicate ratio, 
1.6 and 0.2 respectively. This large difference may not be real, 
even though it seems plausible that the grain growth inside the 
globule is mainly due to gas phase carbon and PAH capture onto 
the graphite grains. 

With the column density of the silicates from the grain 
size distribution, the solar abundance ratio [Si]/[H] {- 3.6-10"^), 
and the assumption, that all silicon are found in the silicates, the 
column densities N(2-//2 -i- H) can be found. For the rim of the 
globule the column density ratio N(2-//2 + H)/Av = 2.3 ■ 10^' 
and towards the more obscured star the ratio is estimated to be 
1.5 ■ 10^'. This is in accord with the average column density 
ratios found for the globule rim earlier (in Paper II). 



4.1.2. The extinction curves 

So far we have only considered the reddening curves, but now 
we can use the models to determine the extinction curves, which 
requires an extrapolation to zero wave number The two models 
give A(K,)/E(I-K,) = 0.34 for star #947 and 0.13 for star #2. 
The difference is striking - but is it reasonable? It is clear that the 
extinction determined towards #2 is more sensitive to systematic 
effects simply because the obscuration is less and thereby small 
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Fig. 10. The excesses from the stars #2 and #947 and their corre- 
sponding model for the whole spectral range. Photometric mea- 
surements are marked with ' H-' signs, while the spectra are drawn 
as full drawn Unes. The models are marked with 'x' signs and red 
lines. 
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Fig. 11. The excesses from the stars #2 and #947 and their 
corresponding model for the wavelength range A > 2 jim. 
Photometric measurements are marked with '4-' signs. The spec- 
tra are drawn as full drawn lines. The models are marked with 
'x' signs and red lines. 



errors due to how well the model atmosphere represents the true 
spectrum of the star, how accurate the photometry is calibrated 
etc will be more important. We also note that there is a small 
decline in the UV in the reddening curve for #2, which does not 
seem very likely (this feature is actually causing the sharp peak 
in the derived graphite grain size distribution). On the other 
hand, it is interesting to note that both these reddening curves 
can accurately be represented by the extinction of only graphite 
and silicate grains with matching size distributions. From this 
point of view, there is no need for another dust component to 
explain the well documented slow decline of the reddening 
curve beyond 2 /^m. The dust model for the rim of the cloud is 
dominated by small grains, which - if it is basically correct - 
means that there is in fact less extinction at and beyond the K 
band than previously thought. Consequently, the reason why the 
reddening curve declines slowly beyond the K band is simply 
because there is not much extinction left. 

As indicated above, the reddening curve for the more ob- 
scured star (#947) is better determined and the corresponding 
extinction model should be trustworthy. In this case, the expla- 
nation why the reddening curve declines slowly beyond the K 
band is the presence of larger (graphite or more general, carbon 
dominated) grains. This is not very surprising. 



4.2. Ices and band strengths 

The integrated optical depths, J Tiv)dv, of the water and CO ice 
bands were calculated. These quantities are directly related to 
the column densities via 

where A is the band strength. iGerakines et al.l (Il995h have in- 
vestigated the band strengths for the stretching modes of O-H in 
water ice at A = 3.02 jum and C = (9 in CO ice at A - 4.67 fim, 
and found that they are not highly dependent on the pollution 
from other substances. The resulting column densities of water 
and CO ice are listed in Table [3] 

4.3. CO gas column density 

Hariunp aa et al.l (|2004|) have observed '^C(9 and C'^(9 emission 
from the B 335 globule and in the neighbourhood of the sightline 
to our background star (A12 in their Table 1) they found a gas 
N(CO)-column density, that deviates from a linear correlation 
with the A j extinction, probably caused by freeze-out of CO. The 
gas column densities were estimated for the two isotopologues 
(Table O, giving slightly different results. 

4.4. Extinction models for ice coated grains 

The ice coatings may either form on the grains by surface re- 
actions or condense onto the surfaces out of the gas phase. In 
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Table 3. Ice column estimates. 



CO column density from submillimeter measurements by Harjunpdd et al (2004) 

ice mass 

isotop N(isotop)10 abundance- 10- @Aj CO gas CO total expected 

[cm^^] column density [pg cm^^] 

Aj = 3.57 as for pos A12 Harjunpdd et al (2004) 

"CO gas ' 45.0 1.108 3.57 19 

''CO total 108.6 1.108 3.57 46 27 



Aj = 3.45 as for sightline towards #947 
'^CO gas 43.5 
"CO total 105.0 


1.108 3.45 
1.108 3.45 


18 

44 26 


Aj = 3.57 as for pos A12 Harjunpdd et al (2004) 
C'**0 gas ' 11.0 0.204 3.57 
C'*0 total 14.5 0.204 3.57 
Aj = 3.45 as for sightline towards #947 

C'*0 gas 10.6 0.204 3.45 
C'*0 total 14.0 0.204 3.45 


25 

33 8 

24 

32 8 


ice column density using laboratory band strength measurements by Gerakines et al( 1995) 

ice bandstrength lO'* jT(v)dv N IO"" ice mass 

[cm] [era"'] [cm^^] column density [/ig cm"^] 


H2O ice 2.00 
CO ice 0.11 
COjice 0.11 


100 61 
3.4 30 
<3 <30 


15 
15 
<20 


ice coat and grain size distribution gained 
from extinction measurements fitted to a modified MRN model 
ice refractive ice mantle volume 10* 
indices [A] [cm] 


ice mass 
column density [pig cm^^] 


H2O ice Fig[13 
CO ice Fig[T5 


Cl 


16 29 
(mixture) 18 21 


28 
23 


silicate band f^; 




column density 


(f5,/volume) 10 ' bandwidth 




volume ■10*' mass 
[cm] [pg cm"^} 


j T(y) dv 
[cm"'] [cm"'] 


#2 0.40 
#947 0.45 




33 115 
44 150 


12 110 
10 110 



" IWarren & Brand! ( l2008h 
' lElsila et al.l ilWh 



Table 4. Column density and mass ratios observed towards star 
#947. 



mass ratio 


column density ratios 




#2 


#947 


ISM" 


graphite/silicate 


0.2* 


1.6' 


0.4 


/f20,c(./silicate 




0.18* 




CO,„/silicate 




0.15* 




COicel HiOice 




1.0* 








0.8'' 




N{CO,„)IN{CO,„,) 




0.8-1 





" Abundance limited ratio (= 1.81) such that all C-atoms exist as 
graphite and all Si-atoms are incorporated into silicates. 
* from modelling 
from band strengths 



both cases it is reasonable to assume that the ice thickness of 
the layers is the same for all grain sizes and thus the ice mass is 
dependent on the size of the grain surface. For w ater ice we use 
the refractive index from I Warren & Brand j (l2008i) and using the 



thickness of the ice layer as a parameter we find a good fit to the 
measurements as can be seen in Fig [13] In these calculations we 
have excluded grains smaller than 25 A, for which the ice man- 
tles evaporate because of the i nstant temperature rise when hit 
by UV photons or cosmic rays dHoUenbach et al.ll2009l) . 

The thickness of the coating in combination with the grain 
size distribution we derived above (see Fig[T2b now allows us to 
calculate the H2O column density. We note that there is more ice 
needed (almost a factor 2, see Table [3]l to explain the absorption 
band when we distribute the ice on the grain surfaces than if we 
simply consider the band strength (which is equivalent to a thin 
sheet of water ice). We have not looked into the detailed reason 
for this difference, but it is probably real and if so, all estimates 
in the literature on column densities of ices are underestimated 
by about a factor two ! 

It is interesting to note that the ice coating slightly changes 
extinction also outside the absorption band and the effect is to 
lower the extinction, see Fig [14] 

For the CO ice the situation is a bit more complicated be- 
cause the dielectric constants are heavily influenced by the 
presence of other frozen molecules, like H2O, O2, N2 as well as 
CO2. Both the peak position and the width of the CO ice band 
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Fig. 12. The grain volume distribution as modelled for the rim 
(probed by star #2) and the central region of B 335 (probed by 
star #947). The grain size distribution model is a power-law dis- 
tribution with simple spherical grains of graphite (full lines) and 
silicates (dashed) with an additional gaussian grain accumula- 
tion at one grain size. 



0.40 - 




0.28 0.30 0.32 0.34 

waveriumber [l//.irr] 

Fig. 13. Detailed view of H2O ice band (fully drawn and a 
model (dashed) of ice c oated grains (with refractive indices from 
IWarren & Brandil2008h . 



Fig. 14. Influence of H2O ice on the extinction in the optical 
wavelength range. The ice (dashed) lowers the extinction in the 
optical from the ice free grain extinction (dotted) hne. 




0.2135 0.2140 

wavenumber [1//.xm] 



Fig. 15. Detailed view of CO ice extinctio n and a model (dashed) 
and the different ices (from lElsila et al.l [T997i) . that have been 
included. 
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are dependent on the mixture and in an attempt to match the 
observed band we have optimised for three different mixtures 
using data from Elsila et al. (1997), see Fig [15] We realize that 
this fit cannot be used to claim that e.g. there must be N2 frozen 
onto the grains - our fit should rather be seen as an example of 
a possible mixture. However, the calculated column density is 
not very dependent on the exact mixture and (like in the case of 
water ice) there is more CO ice needed to explain the absorption 
band when it is located on dust surfaces than achieved from a 
direct calculation using the band strength, see Table |3] 



5. Discussion 

5.1. The reddening curve 

As mentioned, the infrared reddening curve for dark clouds has 
recently been determined in a number of investigations. By as- 
suming the extinction ratio e.g. Ah/Ak to be the same as for the 
diffuse ISM (for which it has been determined), the extinction 
curve follows from: 

Aa/Ak, = (Ah/Ak, - 1)-(Ea-kJEh-k,) + 1- (1) 

As we question (see below) the assumption of a more or 
less universal Ah/Ak value, we have converted the "extinction" 
given in these papers back to the reddening as given in Table 
|5] For comparison we also include the reddening curve for the 
diffuse ISM determined by Indebetouw et al.. (2005). There is 
a good agreement between the different investigations and the 
reddening of our star #2, while our result for star #947 differs 
significantly at the longest wavelengths. For reference we also 
give the reddening produced by the Rv = 5.5 model from lDraind 
(l2003ah . As has become quite clear during the past decade, the 
dip in the wavelength range 3-7 yum, predicted by this model is 
not present in the observations. For this reason, one cannot either 
trust that it can be used as a link between the near-lR reddening 
and the absolute extinction (which is frequently assumed). 

It is interesting to note that the interstellar reddening of the 
diffuse ISM does not significantly differ from that of dark clouds, 
at least until the IRAC 4 band. It is unfortunate, however, that the 
MIPS 24 band has not yet been included in any investigation of 
the diffuse ISM. 

5.2. The extinction 

5.2.1. The diffuse and moderately dense ISM 

In Paper II we studied the extinction of the same dark cloud 
in the optical and near-IR regions based on a large number of 
stars. We came to the conclusion that the CCM functional form 
represents the reddening well. However, this does not neces- 
sarily mean that the extinction curve is adequately described 
by the CCM function at these wavelengths. The problem has 
to do with the extrapolation of the reddening curve to zero 
wave number. Even though variations of the reddening curve 
in the near-IR have been noted, it is generally found that the 
reddening curve can be understood if the extinction in the near- 
IR foll ows a p ower law (A,i oc A^^), where = 1.6 - 1.8 (see 
lDrainell2003bl) . As cited in the Introduction there is, however, 
increasing evidence that such a power law is no longer valid 
beyond 3 yum, and it is therefore not adequate to assume e.g. a 
Ah/Ak ratio based on a power law extrapolation. The obvious 
solution is to extend the observations further into the infrared. 



and this is the approach of the present paper. However, even if 
the Spitzer archive allows us to use observations out to 24 yum, 
the extinction is not zero at this wavelength either. Somehow 
we must extrapolate to zero wave number. For this purpose we 
use a grain model designed to fit the reddening curve all the 
way from the UV to 24/im, including the important silicate 
absorption as well as ice features (for the most obscured star). 
Our grain model is a bit ad hoc (a power law size distribution 
plus a Gaussian distribution centered on a certain size kept 
as a parameter in the fitting process), but it is interesting that 
only silicates and carbon grains (plus ice mantles) are needed 
to explain the reddening curve(s) from the UV all the way to 
24yum, including the region 3-8 yum. Thus, even though we 
cannot claim that our grain models are unique, they enable us 
to extrapolate the reddening curve to zero wave number in a 
physically reasonable way. 

Ideally one would like to measure the relation between 
the absolute extinction and the reddening without being de- 
pendent on a dust model. As discussed in Paper II, star counts 
offer a way to find the true extinction at a given wavelength, 
thereby providing the link between the reddening curve and the 
extinction curve. The main problem in applying this method 
to dark clouds is their small scale structure. In practice large 
uncertainties result. 

Based on a near-IR survey of the Galactic Bulge, 
iNishivama et al.l (l2006l) identify the Red Clump (RC) stars 
in the colour/magnitude diagrams (CMD) for a large number 
of sub-fields. Because of variations of the extinction across the 
field the position in the CMD of the peak of the RC differs 
from one sub-field to another, and the locus of all these RC 
peaks directly gives the A/i-j/E(J-Kj) ratio. Their accuracy is 
very good: Aks/E{J-Ks) = 0.494 +0.006. This is close to our 
result for sta r #2: Ay, /E(J- Kt) - 0.51. It is also interesting 
to note that Fritz et al.l (1201 lb . who determine the extinction to 
the galactic centre by measuring the hydrogen recombination 
lines from 1.28 to 19//m (relating the absolute scale to the radio 
continuum and using Menzel's Case B), get A(Kj)/E(J-Kj) 
«A(2.16)/(A(1.28)-A(2.16)) = 0.51. 

Three different methods thus arrive at the same A(Ks)/E(J- 
K5) ratio. All three investigations probe both the diffuse ISM and 
denser regions, and one may wonder if this r esult is applicable 
for the diffuse and/or dense ISM in general. Ilndebetouw et al.l 
(l2005h investigate two directions close to the Galactic plane (/= 
42 and /= 284) and get a different result: A(Kj)/E(J-K,) = 
0.67+0.07. Their method is also based on the statistical evidence 
that red giant stars have a relatively narrow distribution both 
regarding their near-IR colours and their absolute magnitudes 
and by assuming that the diffuse interstellar matter is rel atively 
smoothly distributed along the lines of sight, Indebetou w et al.l 
(2005) can interpret the colour/magnitude plots in terms of 
extinction ratios and the extinction per kpc. They base their 
analysis on large area surveys, combining 2Mass with Spitzer 
mapping, and the red giant population is well separated from 
the dwarf population in their diagrams (see their Fig. 3). 

Even though the extinction in our line of sight to star #2 
should be dominated by the dark cloud there is also a component 
of diffuse extinction from the path behind the star It is therefore 
of interest to investigate whether e.g. our extinction ratio 
A jI Ak-2.96 would deviate by large f rom the fit to the observe d 
red giant distribution observed by ilndebetouw et al.l (l2005h . 
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Table 5. Reddening (A^-Ax,)/E(J-K,) 



wavelength [fim] 0.35 0.41 


0.51 0.67 


0.80 


1.26 


1.65 


2.16 


3.55 


4.49 


5.73 


7.87 


23.7 


1B335#2 6.91 7.33 
2B335#947 6.18 5.49 


6.91 5.11 
4.73 3.64 


3.82 
2.85 


1.00 
1.00 


0.32 
0.35 


0.00 
0.00 


-0.28 
-0.27 


-0.35 
-0.33 


-0.40 
-0.39 


-0.37 
-0.48 


-0.38 
-0.81 


3 Draine (2003a) 7.15 6.67 
Rv=5.5 


5.74 4.09 


3.07 


1.00 


0.37 


0.00 


-0.38 


-0.49 


-0.56 


-0.55 


-0.60 


4 Indebetouw et al. (2005) 

two large regions close to the Galactic plane 






1.00 


0.36 


0.00 


-0.30 


-0.37 


-0.38 


-0.37 




5 Flaherty et al. (2007) 

star forming regions NGC 2068/2071 






1.00 


0.37 


0.00 


-0.25 


-0.33 


-0.35 


-0.35 


-0.32 


6 Roman-Zuniea etal. (2007) 

dense cloud core 






1.00 


0.40 


0.00 


-0.27 


-0.35 


-0.39 


-0.38 




7 Chaoman etal. (2009) 

molecular cloud (Ophiuchus) 






1.00 


0.37 


0.00 


-0.25 


-0.33 


-0.38 


-0.38 


-0.46 


8 Boosert et al. (2011) 

69 isolated dense cores (from polynomial ibid., equation (1)) 




1.00 


0.35 


0.00 


-0.26 


-0.31 


-0.35 


-0.37 


-0.44 





From their fit they get Aj/Ak = 2.5+0.2 and C/i:=0.15+0.1. 
We use the intrinsic colour index and the absolute magnitude 
of a Kl giant to construct the J versus J-K diagram for such 
stars along a line of sight. By adjustin g the Cy value to 0.105 
(i.e. well within the error estimates of llndebetouw et al]|2005b 
we can exactly reproduce the same curve, see Fig [l6l For the 
corresponding diagram with J versus H-K where we compare 
the curve based on the result of Indebetouw et al. (2005) (with 
Ah/Ak = 1.5+0.1 and Ca:=0. 15+0.1) we get an identical curve 
for our value, AhIAk - 1.63, if we use C/('=0.135. This value 
differs from that we used to make the curves overlap in the 
previous diagram. If we, as a compromise, use the mean value, 
Cy=0.12, we still get cu rves well within the errors quoted by 
llndebetouw et alj ( l2005l) . see Fig.flTl 

What does this comparison tell us? It basically illustrates the 
difficulty to determine the link between the reddening and the 
extinction curves. Not even s uch a careful inves t igation , based on 
high quality data, as that of llndebetouw et al] (l2005l) . excludes 
considerably different results for the diffuse ISM (like our result 
for star #2 or those determined for the path t owards the GC by 
iNishivama etal] (l2006l) and lFritz et all OOllh . 



5.2.2. Dark clouds 

The reddening curve for star #947, sampling dense regions of 
the B 335 globule, differs considerably at the longes t wave - 
length (24^m) from p revi ous work by Flaherty et al.l (l2007l) . 
IChapman etall (l2009l) and 'Boogert et al.' ('201l')(see Table 
On the other hand, the results of these authors agree well to ours 
for star #2, which probes a less dense region than #947. Does it 
also mean that #94 7 traces denser regions than those investigated 
in the cited papers. [Flaherty et al.l (l2007l) analyse data from five 
star formation regions, and there are small differences regard- 
ing the reddening in the IRAC bands. For the MIPS 24yum they 
present data for two regions, NGC 2068/207 1 and the Serpens 
region. Their background sources probe regions with less ex- 




Fig. 16. The black curve represents the locus of typical red giant 
stars at different distances along a line of sight with an extinction 
of Cji-=0.15 kpc~' andE(J-Kj) — 1.5+0.2 A^j as determined by 
llndebeto uw et al.l ([2005') for the diffuse ISM. The dashed curves 
represent the error estimates by these authors. The red curve (on 
top of the full drawn black curve) represents our determination 
of star #2 with c;f =0.105 kpc ' and E(J-Kj) = 1.96 A^,,. 



tinction than that for #947. IChapman et al.l (l2009l) analyse data 
from 3 star formation regions and they look separately at regions 
with different extinctions. Their most obscured sta rs have extinc- 
tions c omparable to that of our star #947. Finallv. iBoogert et al.l 
(201 1"), who analyse a number of dense cores, certainly include 
lines of sight with larger extinction than our case study. They de- 
termine the spectral classes from IR spectroscopy and use stel- 
lar atmosphere models to determine the intrinsic SEDs of the 



12 



S. Olofsson and G. Olofsson: The mid-IR extinction in molecular clouds. 




J-K H-K 

Fig. 17. The same assumptions as in Fig [T7] The dashed curves represent the error estimates by these author s. The red curves 
repres ents our determination of star #2 with Ck-0A2 kpc The red curves are within the error limits quoted bv llndebetouw et aP 
(l200l . This is remarkable in view of the large difference in the E(]-Ks)/Afcs and E(H-K5)/Ajf , betwee n our result for star #2 (wh ich 
probes the outer regions of the B 335 globule plus diffuse clouds behind the globule) and the results of llndebetouw et aT] (l2005h . 



stars. For one of their stars, judged to be well classified, they 
determine the extinction curve for the continuum assuming the 
same AKsfEii-K,) ratio as determined by Indebetouw. In Table 
|5]we have removed that assumption to get the reddening cu rve. 
Con sequently, the reddening curve determined bv feoogert et alJ 
( 1201 1) is only based on one star and their method is the same 
as ours. However, their reddening curve is applied to the 27 
other stars and in most cases it fits well, except frequently at the 
longest wavelength, MIPS 24yum. It is also at this wavelength 
where our results differ significantly. 

There appears to be a problem to include the MIPS 24 fim 
fotometric point in a reliable way in these studies. Not only in 
the l arge sample of individually studied stars by Boog ert et alj 
(1201 1) there are many cases which dif fers significantly from 
the others. Also in the investigation of IChapman et alJ (l2009t) 
the less obscured regions exhibit unrealistic values at 24fj.m. 
Assuming that the technical details (colour correction due to the 
broad filter, the psf correction) have been correctly done in all 
the cited papers, there remains three sources of error: 



1. The MIPS 24fim photometry is not reproducible enough. 
According to the MIPS handbook, the repeatability is 1% 
or better and the absolute calibration is better than 5% 
determined (using observations of solar type stars. A, and K 
giant stars). This cannot be the main reason for the apparent 
problems with the 24 yum photometry. 

2. The model atmospheres do not represent the stellar SED 
well enough for late-type giants. This may well be an 
issue in particular as the model fits are done with relatively 
low spectral resolution, preventing an accurate match of 
metallicity, microturbulence and surface gravity. The most 
important parameter, Tp//, is probably well determined by 
the combination of phot ometry and spe c trosco py. It should 
be noted, however, that iBoogert et al.l (1201 Ih do not use 
high resolution synthetic spectra for their comparison to the 



observed spectra - instead they use opacity sampled spectra 
and degrade the resolution of the observed spectra from R = 
2000 to R = 100 for the comparison. This is in our view a 
waste of information and a likely source of uncertainty. 

3. Late type giants are often irregular variables (and, indeed 
Mira variables). As the observations in the various wave- 
bands are measured at different times, this can cause artificial 
drops and/orrises in the observed SEDs. There is no indica- 
tion in the observations of #947 in B 335 that this is the case 
- actually a smooth SED results even though seven different 
instruments contribute. 

In order to minimise the possible errors due to the stellar 
mo del spectra, w e have plotted the stars in the investigation 
by 'Boogert e t al.l (1201 ih in a diagram with the colour in- 
dex K5-m(24) versus E(J-K.5). The E(J-K5) values are taken 
from their Table 3 in accordance with their assumption, E(J- 
Kj)=1.5 A(K5). In the same graph (Fig.fTSl) we also include our 
star #947, and it is clear that the position of this star in the plot 
significantly differs from the general trend, which is actually 
well determined. If we neglect the (small) difference in the 
intrinsic colour index, (K5-m(24))o, for the individual stars we 
get 

E(K,-m(24)) = (0.357+0.013)E(J-K,) 

This value is lower than iBoogert et al.l (1201 ih derived from 
one star in their sample (0.44) and much lower than our value 
for star #947 (0.84). 

Is there a reasonable explanation for this difference? We 
first note that the estimated error due to classification uncer- 
tainty and the photometric accuracy is within the size of the 
symbol of #947 in Fig. [18] so it should be excluded to explain 
the exceptional position of this star in the diagram by poor 
photometry or classification. Even though B 335 is not very far 
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Fig. 18. This c olour/colour diagram is based on data from an 
investigation bv' Boogert et alJ (1201 ih of field stars behind cloud 
cores. The slope of the curve is well defined, 0.357±0.013. The 
asterisk represents the star #947 in our investigation, and it is 
clearly different from the general trend. 

fro m the galactic plane (b = -6.5), most of the cores investigated 
bv lBoogert et alJ OOllI) are located closer to the galactic plane 
and/or towards the galactic center. This means that much of 
the continuum extinction can be due to the path behind the 
cores and as a consequence much of the determined E(J-K5) 
could be due the diffuse ISM, while most of the mid-IR excess, 
E(K5-m(24)) could be caused by the dense cores. Such an 
explanation might be plausible, but we cannot rule out that our 
case study of the extinction towards star #947 behind a dense 
part of B 335 reflects unusual cloud properties. 

6. Conclusions 

Based on photometry in the range 0.35 - 24 jum as well as IR 
spectroscopy we determine the reddening curve towards stars 
behind the dark globule B335. For the most obscured star we 
also determine the absorption bands due to water ice, CO ice 
and sihcates. - We confirm earher findings that the reddening 
curve levels ofi^ beyond the K band in contrast to current models 
(e.g. IWeingartner & Draind 1200 1). From the reddening curve 
we also draw the following conclusions concerning the grain 
size distribution, the extinction and the ice bands. 

Grain size distribution. 

The grain size distribution features. Based on a relatively simple 
assumption for the grain size distribution (a power law and 
a Gaussian) for silicate and graphite grains we find solutions 
that closely match the observed reddening curves for the two 
background stars for which we have silicate band observations. 
One of the stars probe the central part of the globule and the 
other the rim. For the most obscured star the model fit has a 
peak of graphite grains at a radius of 2 fim, while the graphite 
peak for the other star is at 0.5 fim. In neither case there is any 
need for an additional grain material to explain the reddening 
curve. 



The abundance ratio graphite/slicates differs substantially 
for the two models (much more graphite in the central region). 
As a consequence we have two different explanations of the 
shallow reddening curve beyond the K band: for the denser 
region it is because of larger graphite or, more general, carbon 
grains which give rise to a higher mid-IR extinction than 
predicted by previous models constructed to match the UV to 
near-IR spectral range (cf. Weingartner & Draine 2001). For the 
outer region of the globule the explanation is simply that the 
extinction as such is low beyond the K band. Actually, much 
of the extinction towards the stars at the rim of the globule 
may originate in the diff'use medium beyond the cloud, and we 
suggest that this explanation to the shallow reddening curve 
beyond the K band is valid not only for the outer regions of the 
cloud but also for the diffiise ISM. 

Gas column density. The column densities ratios 
Nil ■ H2 + H)/Av are estimated to be 2.3 • 10^' for the 
lower extinction sightline and 1.5 • 10^' for the higher extinction 
sightline. 

Extinction model. 

The extinction model for the outer part of B 335 gives the 
relation Aks = 0.51 ■E(J-Ks). This latter relation agrees 
with recent determinations towards the galactic centre by 
Nishivama et al. (2006) and Fritz et al. (2011) using different 
methods. ilndebetouw et al.l (l2005i) get a higher value for the dif- 
fuse ISM, Aks - 0.67-E(J-Kj), but we show that this diff'erence 
is not significant and conclude that A^j x 0.51 •E(J-K5) for the 
diffuse ISM as well as for moderately dense molecular regions. 

The extinction model for the denser part of the B 335 
gives the relation Aks - 0.97 ■E(J-Kj). Although this model 
describes the reddening curve from the UV to 24 yum, including 
the silicate, the water-ice, and the CO-ice bands well, we have 
indications that it represents conditions that may be different 
from other dense cores. 

Ice bands. 

The water ice band is reproduced by grains coated with pure 
water ice. The column d ensity of wate r ice is half of that found 
in the Taurus clouds (IWhittet et al.l (l200lh 's) for the same 
extinction. The difference is likely due to a larger interstellar 
radiation field in the B 335 region, desorbing ice mantles deeper 
in the cloud. 

The CO ice band differs in position from that of pure CO 
ice mantles, indicating that CO is mixed in other ices. Based on 
published laboratory measurements of frozen gas mixtures we 
show that the CO ice band can be understood if CO is part of 
such mixtures (where H2O is a natural part and N2 a possible 
part). The ice/gas ratio is relatively high (0.8 - 1.) and as this 
is an average along the Ay = 3.5 cloud extinction towards the 
probing star and as we must allow for relatively deep regions 
free from CO ice on both the front and the rear side of the cloud, 
the ice/gas CO ratio must be greater than unity in the interior 
of the cloud, and more t han half of the gas i s depleted. This 
confirms earlier findings (Bacmann et al. 2002) that CO (or CO 
isotopologues) cannot be used as a mass tracer in dense and 
cold molecular clouds. 

The CO2 ice band at 15.2 p.m indicated in the ISOCAM- 
CVF spectrum is not confirmed in the Spitzer spectra and we 
can only give an upper limit to the COa-ice column density. 



14 



S. Olofsson and G. Olofsson: The mid-IR extinction in molecular clouds. 



Acknowledgements. . 

This publication makes use of data products from the Two Micron All 
Sky Survey, which is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California Institute of Technology, 
funded by the National Aeronautics and Space Administration and the National 
Science Foundation. 

This work is based [in part] on archival data obtained with the S pitzer 
Space Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. 



References 

Bacmann, A., Lefloch, B., Ceccarelli, C, et al. 2002, A&A, 389, L6 
Boogert, A. C. A., Huard, T. L., Cook, A. M., et al. 201 1, ApJ, 729, 92 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Chapman, N. L., Mundy, L. G., Lai, S., & Evans, N. J. 2009, ApJ, 690, 496 
Chiar, J. E., Pendleton, Y. J., Allamandola, L. J., et al. 201 1, ApJ, 731, 9 
Draine, B. T. 2003a, ArXiv Astrophysics e-prints 
Draine, B. T. 2003b, ARA&A, 41, 241 
Draine, B. T. 2003c, ApJ, 598, 1017 
Draine, B. T. & Lee, H. M. 1984, ApJ, 285, 89 
Elsila, J., Allamandola, L. J., & Sandford, S. A. 1997, ApJ, 479, 818 
Flaherty, K. M., Pipher, J. L., Megeath, S. T, et al. 2007, ApJ, 663, 1069 
Fritz, T. K., Gillessen, S., Dodds-Eden, K., et al. 2011, ApJ, 737, 73 
Gerakines, P. A., Schutte, W. A., Greenberg, J. M., & van Dishoeck, E. F. 1995, 
A&A, 296, 810 

Gibb, E. L., Whittet, D. C. B., Boogert, A. C. A., & Tielens, A. G. G. M. 2004, 
ApJS, 151, 35 

Harjunpaa, P, Lehtinen, K., & Haikala, L. K. 2004, A&A, 421, 1087 
Hauschildt, P H., Allard, F, Ferguson, J., Baron, E., & Alexander, D. R. 1999, 
ApJ, 525, 871 

HoUenbach, D., Kaufman, M. J., Bergin, E. A., & Mehiick, G. J. 2009, ApJ, 690, 
1497 

Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 931 
Knez, C, Boogert, A. C. A., Pontoppidan, K. M., et al. 2005, ApJ, 635, L145 
Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Nishiyama, S., Nagata, T, Kusakabe, N., et al. 2006, ApJ, 638, 839 
Olofsson, S. & Olofsson, G. 2009, A&A, 498, 455 
Olofsson, S. & Olofsson, G. 2010, A&A, 522, A84+ 
Pontoppidan, K. M., Fraser, H. J., Dartois, E., et al. 2003, A&A, 408, 981 
Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618 

Roman-Zuniga, C. G., Lada, C. J., Muench, A., & Alves, J. F. 2007, ApJ, 664, 

357 

Skrutskie, M. P., Cutn, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 

Warren, S. G. & Brandt, R. E. 2008, Journal of Geophysical Research 

(Atmospheres), 113, 14220 
Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 

Whittet, D. C. B., Gerakines, P. A., Hough, J. H., & Shenoy, S. S. 2001, ApJ, 
547, 872 



